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We present a general framework for contextuality tests in phase space using displacement oper¬ 
ators. First, we derive a general condition that a single-mode displacement operator should fulfill 
in order to construct Peres-Mermin square and similar scenarios. This approach offers a straight¬ 
forward scheme for experimental implementations of the tests via modular variable measurements. 

In addition to the continuous variable case, our condition can also be applied to finite-dimensional 
systems in discrete phase space, using Heisenberg-Weyl operators. This approach, therefore, offers 
a unified picture of contextuality with a geometric flavor. 

PACS numbers: 07 . 10 .Cm, 03 . 65 .Ta, 03 . 65 .Ud 


The concept of contextuality highlights the fundamen¬ 
tal departure of the quantum description of the world 
from that of classical models. In everyday life, the re¬ 
sult of a measurement occurs irrespectively of which 
other compatible measurements (the so-called “context”) 
are simultaneously measured together with it. Quan¬ 
tum mechanics, however, as originally shown by Kochen 
and Specker Dia, does not allow the noncontextual as¬ 
signment of values to physical quantities, giving rise to 
testable differences in predictions of quantum and non¬ 
contextual theories [5], as in Bell’s theorem [1] (locality 
being a special case of noncontextuality). Remarkably, 
this can happen irrespectively of the quantum state, i.e., 
state-independent contextuality (SlOfKUIO]. Contextu¬ 
ality has been confirmed in experiments with genuine mi¬ 
croscopic systems, such as two qubits and qutrits 

[HHsi, and it has been recognized as a critical resource 
for quantum computing [Mm and for several quantum 
information-theoretic tasks [2(11122) . Recently, there has 
been a growing interest in exploring quantum features 
such as quantum superposition closer to the classical 
realm. Along this line, a far-reaching attempt, perhaps, 
is to experimentally probe the quantum contextuality in 
continuous variable systems with most “classical-like” op¬ 
erations in phase space. 

Phase space displacement operators exhibit favorable 
properties, making them very suitable for investigating 
contextuality in phase space. For one thing, their phase 
space functions given by Weyl-Wigner correspondence 
are bounded and, thus comply with the explicit assump¬ 
tion used for the derivation of the existing noncontextu¬ 
ality inequalities in which the classical variables take on a 
limited range of values. This is in contrast to, e.g., parity 
operators as they are described by an unbounded delta 
function in phase space. That is why parity correlation 
measurements can lead to the violation of Bell’s inequal¬ 
ity for states represented by positive Wigner functions 
[23] ■ Moreover, displacement operators, unlike quantum 
operations such as Kerr nonlinearities, do not transform 
coherent states into nonclassical states characterized by 


negative-valued Wigner functions. Indeed, displacement 
operators reveal a characteristic geometrical feature of 
quantum mechanics in phase space originated from the 
canonical commutation relations between position and 
momentum operators. 

In this Letter, we present a general approach to quan¬ 
tum contextuality in phase space, and in particular to 
SIC. This is done by giving an explicit recipe to wit¬ 
ness SIC in continuous variable systems. Our results in¬ 
clude as a special case a previous proposal of observing 
quantum contextuality in infinite dimensions [24) . Our 
approach offers a simple and optimal scheme for experi¬ 
mental implementations of the tests via measurements of 
modular variables, i.e., the Hermitian components of the 
displacement operator. Modular variables play a signifi¬ 
cant role in various tests of quantum nonclassical effects, 
such as the GHZ theorem [5^, macroscopic realism [2S] 
and entanglement detection mm- 

Preliminary notions. — Let us consider nine di¬ 
chotomic observables Aj^, j, fc = 1,2,3, such that triplets 
of observables sharing a common index j or k are com¬ 
patible. Here, compatibility means that they can be 
measured jointly or sequentially without disturbing each 
other. In quantum mechanics this is, for instance, the 
case if the observables commute. Thus, we can measure 
the following mean values [3: 

(XPm) =(^ll7ll2^13) + (^21^22^23) + (^31^32^33) 

+ (A11A21A31) J- (A12A22A32) — (A13A23A33). 

( 1 ) 

Such conditions are satisfied by the following choice of 
{Ajfe}, also known as Peres-Mermin (PM) square [^130) . 

Ajk k = 1 k = 2 fc = 3 

j = 1 cTj, (g) 1 az®iJz ^2^ 

j = 2 l^CTx Ua; 0 1 CTa, 0 (Ta, 

j = 3 (Tz 0 (Ta; CTa: 0 Uz <7y 0 CTy 

leading to a value of (xpm) = 6, for any quantum state, 
since the product of operators along each row and column 
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FIG. 1: Pictorial representation of three displacement ampli¬ 
tudes Oi which span a triangle of area ^ = 7r/4 and therefore 
satisfy the condition in Eq. § used for constructing a PM 
square in phase space. 


is 1, except in the last column, where it is —1. On the 
other hand, the maximal value of (xpm) is 4 for all non- 
contextual hidden variable (NCHV) theories, i.e., classi¬ 
cal probability theories assigning a definite ±1 value to 
each observable independently of the measurement con¬ 
text. This can be verified by a direct substitution of all 
possible ±1 values. Hence, the PM square exhibits SIC. 

To obtain the same construction in phase space, we in¬ 
troduce the single-mode displacement operator, 'D{a) = 
goat-a where (a) is the creation (annihilation) op¬ 
erator of a single bosonic mode and a is a complex dis¬ 
placement amplitude. By using the Baker-Campbell- 
Hausdorff formula, i.e., where 

[X, Y] is a constant, we obtain the known relations 

'D{ai)'D{aj) = -I- a^), (3) 

and 

V{ai)V{a,) = )!?(«,)• (4) 

The phase factor indicates that displacements in different 
directions do not commute in general. 

Phase-space PM Square. — Geometrically, it is illumi¬ 
nating to think of the amplitudes a as real vectors in 
a two-dimensional space, i.e., a := (Rea, Ima), which 
gives 

Imjaia*} = a^ x Uj := |ai||aj| sin6>y, (5) 

with Oij the angle between the two vectors. Now, con¬ 
sider a set of amplitudes ai,a 2 , and as satisfying the 
constraint 

ai X a2 = 02 X as = as X ai = ±^, (6) 

from which ai -|- 02 + 03 = 0 follows. This constraint 
yields the following algebraic relations for the displace¬ 
ment operators similar to Pauli matrices, namely, 

V{ai)V{a2) = ±iV{-a^) , {V{ai),V{a,)} = 0. (7) 


A simple example satisfying condition ( [^ is g iven by am¬ 
plitudes Oi with equal lengths \ai \ = ~ 1.34 and 

rotated by an angle 27r/3 with respect to each other. As 
we will show below, this condition also extends to the 
discrete phase space. 

From Eq. Q we can straightforwardly construct a PM 
square analogous to © for two bosonic modes, 

Ajk _ k = l _ k = 2 _ fc = 3 

j = l Vi{-ai) T>2(-ai) I?i(ai)T>2(ai) 

j = 2 T>2{—0'2) I?i(— a2) T>i(a2)T’2(a2) 

j = 3 I?i(ai)I?2(a2) I?i(a2)T’2(ai) 'Di{az)V2{oLz) 

( 8 ) 

where I?i(I? 2 ) denotes a displacement operator for mode 
1(2). Operators Ajk = -|- within each row or 

column are mutually commuting, and the same holds for 
their real and imaginary Hermitian parts, i.e., the modu¬ 
lar variables and The product of three operators 
in each row and column is 1, except in the last column 
where it is —1. The minus sign arises as a geometric 
phase proportional to the area A = 7 r /4 covered by the 
corresponding loops in phase space, cf. Fig. 

The displacement operators in (|^ do not represent 
physical observables, and in each experimental run, we 
can only measure either the real or imaginary parts, Aj^ 
and A^j,. Therefore, instead of Eq. Q, we now consider 
the real-valued form 


(xpm) = {Ri) + (i?2) + (R3) + (Cl) + (C2) - (C3), ( 9 ) 

where Rj = Re(AjiAj 2 Aj 3 ) for products of three dis¬ 
placement operators in the same row and, analogously, 
Cfc = Re(AifcA 2 feA 3 fe) for operator products in the same 
column. Explicitly, in terms of modular operators we 
obtain 

n. _ ( aR aR _ aI aI \aR _ ( aI aR I aR aI \ aI 
r<^ — ( aR aR _ aI aI \ aR _ ( aI aR 1 aR aI ^ aI 

— [^lk^2k ^lk^2k)^3k \^lk^2k + ^lk^2kl^3k- 

( 10 ) 


Classical bound .— Since only A^. or Aj^ can be mea¬ 
sured, and not directly Ajk, we must compare the ob¬ 
served values with the most general NCHV theory, where 
those observables are associated with independent vari¬ 
ables. In this case, the products Rj and Ck can assume 
values between —2 and - 1 - 2 , and for correlations measured 
as in Eqs. (10), we obtain an upper bound of 12, which 
is the double of the quantum value (xpm)qm = 6. To re¬ 
cover a lower classical bound, we enforce the additional 
constraints ((Aj|)^-b (Aj^)^) < 1, satisfied in QM, to our 
NCHV model, by adding “punishment” terms that lower 
the classical bound whenever such a condition is violated 
(see, e.g., [a 1311132]). Thus, we can define 


XpM — Xpm — + i^jk)"^ ~ l| (11) 

jk 
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and proceed to maximize the expectation value of Xp^ 
with the only assumption that and are rep¬ 
resented by noncontextual classical random variables 
taking values in [—1,1]. It can be proven (cf. 
Appendix) that for A > 2, the maximal value for the 
expression 0 in any NCHV theory is 3-\/3- This ex¬ 
tends the bound derived by Plastino and Cabello 
which holds only for a restricted set of NCHV theories 
where the condition {{A^f^Y + < 1 is assumed, 

and makes contextuality tests in phase space generally 
applicable. 

Contextuality with a single mode. — In close anal¬ 
ogy to the case of qubits, the PM square Q has been 
constructed for two modes. A natural question arises, 
whether a similar set of operators with identical rela¬ 
tions could already be identified for a single continuous- 
variable degree of freedom. As we show now, this is not 
possible in a strict sense, but a PM square for single 
mode displacement operators can be realized with an ar¬ 
bitrary good approximation, when large displacements 
are allowed. 

Let us consider a set of nine single mode displacement 
operators Ajk = T>{Ajk), labeled by the correspond¬ 
ing two-dimensional displacement vectors Ajk- To ob¬ 
tain the same relations as for the two mode case, we 
require first of all, that the product of three operators 
within each row and column commute and are propor¬ 
tional to 1. This is given by the conditions Aji x Aj 2 = 
Aj2 X Aj3 = Aj3 X Aji = kn, etc, for some integer k, 
implying Aji + Aj 2 + Aj^ = 0, etc. Finally, the product 
of all three displacements is - 1-1 meaning that k must be 
even, except for the third column where the product is 
— 1, and therefore, k must be odd. By putting together 
the above relations for the displacement vectors we prove 
the following result: 

Observation. It is not possible to represent the com¬ 
mutation and product relations of the PM square with 
displacement operators on a single mode. 

By the linear relations, it is sufficient to consider four 
vectors An, A 12 , A 21 , A 22 , and define A 13 = -An - A 12 
etc. Then, the cross product relations give us 

All X A12 = kiTT, A21 X A22 = k2'!T^ A31 X A32 = k^TT, 
All X A21 = k^TT, Ai 2 X A22 = Astt, Ai3 X A23 = kgTT, 

where fci, ...,^5 are even, while fcg is odd. A direct, but 
tedious, calculation allows us to remove the vector vari¬ 
ables, obtaining a conditions only for the integers ki 

iki - k 2 f + {ki - knf -f {k 2 - knf - kj - k^ - kj = 
(fc4 - k^f + {ki - kef + (fcs - kef -kl-kl-kl 

From this, one can directly see that, if the first five ki 
are even, the remaining ke has to be even. The other 
interesting cases, where, e.g. three ki are even and the 



FIG. 2: Illustration of the displacement amplitudes Ajk used 
for the approximate construction of the PM square for a single 
mode. See the text for further details. 

three remaining kj are odd, are also not possible. This 
proves the claim. 

Nevertheless, as we will show now, it is possible to ap¬ 
proximate the PM square on a single mode arbitrarily 
well, if the displacement is allowed to be large. Since any 
experimental implementation of the perfect PM square is 
only an approximation, this can still be a way to test hid¬ 
den variable models, if the occurring errors are quantified 
by additional measurements of the compatibility |31j . 

To start, we set the vectors of the first row as 



implying that the length of A 13 is L = f2l. The other 
rows are defined in a similar way, but rotated by 27r/3 (cf. 
Fig.[|). Again, we have to impose the conditions on the 
area of the spanned triangles, but here we make a slightly 
different choice, by demanding that the product of five 
rows or columns should be — 1 , while for the remaining 
one it should be -1-1. One can easily see that this leads 
to the same contradiction as the usual PM square. 

Our conditions can then be written as follows: 

TT 

with odd ki for i = 1,2,3, corresponding to the three 
rows. For the first and second column, it follows that 



with odd ki for i = 4,5. The factor a/3/4 comes 
from the fact that the vectors form an equilateral tri¬ 
angle (see Fig. [^. Finally, for the last column, we have 

\l % = ke^ with an even ke- This, however, follows au¬ 
tomatically from the previous condition, since = 2 £^. 
With odd r and s big enough, we can approximate arbi¬ 
trarily well a/3/4, i.e.. 
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We then choose (. ss I = v^, and the conditions in 


Eq. (131 are automatically fulfilled. Furthermore, with 


this choice we have 

e_ 

2 


3 TT 

4 2 


(16) 


so the conditions from Eq. (14) are approximately ful 


filled. The approximation improves with increasing val¬ 
ues of r and s, which requires an increasing displacement 
i. To give a concrete example, with displacements up to 
a maximal |a| = 10, the condition is L = < 10 al¬ 

lows for a value of s = 15. Then, the choice r = 13 leads 


to an approximation of Eq. (141 with a value ki = 12.990 


instead of the ideal ki = 13. This is already a very good 
approximation of the PM square. 

Discrete phase space. — The generality of our approach 
also allows us to investigate contextuality also in dis¬ 
crete phase space, i.e., a particle taking on discrete po¬ 
sition values with a finite number of position sites. The 
single-step momentum and position displacements in a 
d-dimensional Hilbert space %, respectively, are [33], 


^ _ ^i^TzQld ^ _ ^—i2'n:Pjd 


(17) 


with Q = P = the dis¬ 

crete position and momentum operators, respectively. 
The operators ( [T7| ) are known as Heisenberg-Weyl (HW), 
or generalized Pauli operators. The position and momen¬ 
tum basis states are related via discrete Fourier transform 




|n), (18) 


where the sum is defined modulo d, and they obey 
X"^\n) = \n + m) and Z'‘\k) = \k + 1). The periodic 
boundary condition gives Z‘^ = 1 and = 1. The 
commutation relation between Z and X operators are 


^i27rlm/d 


(19) 


and a general displacement operator can be defined as 

V{l,m) := (20) 

describing the displacement of position and momentum 
with m and I steps, respectively. With the shorthand 
notation aj = i/Tr/d{lj — irrij), then, Eqs. (19) and (20) 
give a product rule for HW operators analogous to Eqs. 
(§ and 0 - The condition established in the continuous 
limit, i.e., Eq. 0 , for the discrete scenario then becomes 


7712^1 — TO 1/2 = rnzh — W 2/3 = = d/2. (21) 


Since the left hand side of this equation is always an in¬ 
teger, the dimension d must be an even number. As an 
example, for d = 2 , the hopping steps are mi = 0 , Zi = 
l,m 2 = 1, /2 = 0, m 3 = 1, ^3 = 1 giving 'D{ai) = cr^. 



FIG. 3: A quantum circuit composed of a sequence of three 
Ramsey measurements used for measuring the correlations 
between modular variables. 


D{a 2 ) = Ox and D{a^) = —iOzOx = Oy, respectively, and 
thus the PM square (|^ for two-qubit systems is recov¬ 
ered. 

Implementation .— The experimental violation of the 
classical bound for requires the correlation mea¬ 

surements of modular operators appearing in (10). By 
adapting techniques described in Refs. m iMi, this 
can be achieved in current trapped ion experiments and, 
eventually with macroscopic mechanical resonators also. 
In general, we consider a qubit-oscillator model whose 
Hamiltonian reads (?i = 1) 


H = cco|t)(tl + + X{tp + \*{t)b |t)(tl, (22) 


where Wg is the qubit splitting and the \{t) are time- 
dependent couplings which, in trapped ions, for example, 
can be realized by state-dependent optical forces (see, 
e.g., |35l|36]). The mode b is the linear combination of 
two (spatial) orthogonal modes: b = cos(d)ai -|-sin( 0 )a 2 . 

Now, we consider the sequence of three Ramsey mea¬ 
surements shown in Fig. ([^. At the beginning of each 
sequence i = 1,2,3, the qubit is first prepared in a su¬ 
perposition (I I) -I- t))/'s/2 by applying a fast 7r/2 

rotation R^{(j)i). In a frame rotating with the bare qubit 
frequency ujq the coupled qubit-oscillator system then 
evolves under the action of H for a time Tt 


U{n) = [|mi + e*^("-)l?b[a(r.)]|t)(tll C/o(r), (23) 


with C/o('r) = being the oscillator’s free evolu¬ 

tion. In addition, the system accumulates phase <(>(rj). 
Finally, the qubit is rotated back by applying i?.n./ 2 ( 0 ), 
the population difference Z is measured and the qubit is 
reset into the state | |). As shown in more detail in [26j . 
the average outcome of this measurement is 

(ZtJ = {Q{ipi,ati)} = (cos {ipi + atp + a^.b)), (24) 

where aj. = a(Ti)e®‘^‘‘ and ipi = (j)i + 4>{Ti). Therefore, 
by choosing either = 0 or = 7 r /2 and by appropri¬ 
ately adjusting ai and 6 (for example, by changing the 
laser direction), we can measure (Aj).) and (Ajj,). Fur¬ 
ther, in the case of interest where different displacement 
operators commute, it can be shown [26| that 

{Zt^Zt^Zt^) = Tr{Q((pi, at J(3((p2,at2)Q(7’3,at3)po}, 

(25) 
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which then allows one to measure the operator products 


Rk and Cj given in Eqs. (101 in terms of correlations 


between the outcomes of three Ramsey measurements. 
In summary all terms appearing in the expression for 
Xp™ can be measured in this way. The symmetric choice 
of single-mode amplitudes given below Eq. Q requires 
displacements of only \ai\ ~ 1.34, which is well within 
what is achievable in current ion trap experiments. As 
an explicit example we provide a detailed description of 
the implementation of C3 in the appendix. 

Conclusions and outlook .— The geometric feature of 
quantum mechanics demonstrated by phase space dis¬ 
placements leads to the impossibility of a noncontextual 
assignment of measurement results. We derived a gen¬ 
eral condition, which is sufficient for constructing vari¬ 
ous SIC scenarios both in continuous as well as in dis¬ 
crete phase space. We showed that displacement (or 
Heisenberg-Weyl) operators offer a unified geometric pic¬ 
ture of contextuality. Moreover, our approach presents a 
simple and symmetric choice of modular variables, sim¬ 
plifying, considerably, the experimental implementations 
using the state-of-art techniques of trapped ions. 

In light of this, it would be interesting to develop our 
approach to Bell nonlocality tests, and construct Bell-like 
inequalities for probing nonlocality with suitable single¬ 
mode displacement operators. Einally, it can be shown 
that the geometric phase A giving rise to contextuality 
diminishes under decoherence and disappears in classi¬ 
cal limit. This motivates us to consider it as a quantum 
resource which can be used for characterizing and quan¬ 
tifying contextuality. 
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APPENDIX 


The appendix provides detailed derivations of the classi¬ 
cal bound and the physical implementation of the mea¬ 
surements involved in the tests. 


(a;i,...,a;9 ,yi,...,?/ 9). Eq. 0 can then be rewritten 
as 


X^pZi^)=FiX)-XPiX), (26) 


where E is a function giving the PM expression of 
Eq. (§ and P the function giving the punishment 
term. The classical bound is given by max^ where 

S = {I^Cil < 1, \yi\ < 1}. We then split our domain as 
S — B UTZ, where 

B = {|x? + y^\ < 1 }, n = {\x^, + yl\ > 1} n 5. (27) 


We can easily compute maxg x?w follows. Since E is a 
harmonic function, i.e., AE = 0, and E is a compact set, 
i.e. closed and bounded, the maximum of E is achieved 
on the border dB. On the other hand, the minimum of 
P is clearly achieved on dB and it is 0. Such a value for 
E has been computed in Ref. [24]. 

To complete our proof, we need to show that for ap¬ 
propriate values of A, this is also the maximum of XpM 
on S. First, we notice that we can get rid of the ab¬ 
solute value in E, i.e. E(X)|.,j = HAjp — 9. We can 
then compute the partial derivative of x^mO^) E^long 
the direction X^ = (0 ,..., 0, cci, 0 ,..., 0 , 0 ,..., 0)., i.e. 

Vx^pmO^) ■ A simple calculation shows that 


VxrM(X).X,/|lX,;|| 


Rj + C^ x^i+yf 

yf V xf + yf ’ 


(28) 


which is always non-positive for 


A > 


max(i?i + Ci) 
2min(a;f -f yf) 



(29) 


Notice that Ri + Ci, as well as E, can be maximized 
on the cube S by simply taking the maximum over all 
values on the vertices of the cube, since it is a harmonic 
function. 

The last thing to prove is that every point in the cube 
can be reached from a point in B with finite increments of 
/i^Xp i.e., by moving in directions where partial deriva¬ 
tive is non-positive. This follows from the fact that for all 
X € 5, there exist Xq S B and numbers fXi > 0 and vec¬ 
tors Xi = (0,..., 0, Xi, 0,..., 0, yi, 0,..., 0), i = 1,..., 9, 
such that 


X = Xo-f ;^/r,X,. (30) 

i 


I. EXPLICIT COMPUTATION OF THE 
CLASSICAL BOUND 

To simplify the notation, let us denote the clas¬ 
sical variable associated with and Aj/,, respec¬ 
tively, as X 3 k+j and ysk+j, and define the vector X = 


II. EXPLICIT DERIVATION OF THE 
THREE-POINT CORRELATIONS 


Here, we provide a detailed derivation of Eq. (25) fol¬ 
lowing the general scheme presented in Ref. [2S] . At the 
beginning of the Ramsey sequence the qubit is initialized 
in state | j,). The pulse sequence Urm{t) acts on the 
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composite qubit-resonator system. Here the Rtt/2{4>) de¬ 
note 7r/2-rotations of the qubit with an adjustable phase 
(j), which in the basis {| j,), | f)} is dehned as 





(31) 


The evolution between the pulses, C/(t) 


I4')('H + 


e®'^('^)l?;,(a(r))|t)(t| Uo{t), describes the qubit-resonator 


interaction for an interaction period r. If at the initial 
time to the resonator is in an arbitrary state po, the to¬ 
tal system density operator after the pulse sequence is 
^i)I4 ')(4'I ® PoUl,[{p, ti) and after partial trace over 
the qubit the oscillator’s reduced density operator reads 


Pih) =P+P^{ti) +P-P {ti). (32) 


The probabilities p-|_ and p- for finding the qubit in state 
I t) and I I), respectively, are then given by 


^±(^ 1 ) =T^T^{E±{p,n)PoEl{ip,Ti)}. (33) 

where E±{ip,T) = ^ [l ± e*‘^('’’^I?(a(r))] Uo{t) are Kraus 
operators satisfying E^_^_E^ + e'^_E- = 1. Right after the 
sequence a projective measurement of Z is done on the 
qubit with expectation value given by 

{Zt^) = p+{ti) - p-{ti) = TT{E+{ipi,Ti)poEl{ipi,Ti)} 

- TT{E_{ipi,Ti)poEi{ipi,Ti)} = Tr{Q((/3i,atJpo}- 

(34) 

Depending on the measurement outcome the conditioned 
resonator state is 

±,, , E±{ipi,Ti)poEl{ipi,Ti) 

P (ti) = -TTX-• (35) 

P±[ti) 

We now generalize the above considerations for a se¬ 
quence of measurements. We set to = 0 the time right be¬ 
fore the first measurement and denote by tn the time after 
the n-th RM is complete. The variables Z{tn) = ±1 de¬ 
scribe the outcome of the respective measurements. Each 
Ramsey sequence is characterized by displacement am¬ 
plitudes a(T„), an adjustable phase of the first 7r/2-pulse 
ipn and the geometric phases (j)n = ipniTn)- Starting from 
the initial resonator density operator po = p(0) the state 
conditioned on the first measurement outcome 771 = ± is 

p''i(tl) = Er„iTi)poEl^{Ti)/p^^ = 

^;^{po{h)+'D{o!i)po{ti)'D{-ai) ± ['D{ai)po{ti) + H.c.]). 

(36) 

This state evolves freely for a time t 2 — T 2 — ti and then 
a second measurement is performed. By repeating the 
arguments from above, the probabilities for this second 


measurement, conditioned on the first outcome, are given 

by 


P±\7]l — 

^ ± ^Tr{Q((752,a(T2))Do(t2 - ti)p^^ {ti)ul{t2 - h)}. 

(37) 

Therefore, the conditioned expectation value of the sec¬ 
ond measurement is 

{Et 2 )vi =P+\vi-P-\vi = Tr{Q((^ 2 ,a(r 2 ))p’'Hi 2 )}, (38) 

where p^^{t 2 ) = Uo{t 2 — ti)p\^ {ti)Ul{t 2 — ti) denotes 
the time evolved conditioned density operator. The two 
point correlation function between two successive mea¬ 
surements is 

{Zt^Zt^) = {p+\+ - p-\+)p+ + (p_|_ -p+\-)p- 
= P+{Q{P2 ,oi{t 2)))+ - P-{Q{p2,a{T2)))- 
= Tr{(5((p2,a(T2))[t/o(t2 - ti)e^'f‘^'D{a{Ti))p{ti)ul{t2 - h) 
+ H.c.]} 

= ^r{Q{Pi,atJQ{(p2,at2)po}, 

where we assumed [T>{atj^),'D{at 2 )] = 0. For a 3-point 
correlation function we have 


{ZtsZt^Zti) — p+{ti){Zt2Zt2)+ — p-{ti){Zt2Zt2)-- 


(39) 

Simply using the above argument we have 
(^t 3 ^t 2 )± = Ev{Q{^p2,at2)Q{p3,at2)p^{0)}, where 

p^{Q) = UQ{ti)p^{ti)Uo{ti). Therefore, we have 

(Zt^Zt^Zt^) 

= o:ti)Q{p2, at^)[Ul{ti)e'-‘^^V{ai)poiti)Uoiti) 

+ H.C.]} 

= xTr{(3((^i, at JQ((p 2 , ai2)[e*‘^^25(ati)Po + H.c.]}. 

(40) 


For mutually commuting displacements ^(afj) 

and Vlat^) then, we obtain 


(Zt^Zt^Zti) = ^^:{QiPl,at^)QiP2, OitJQiP2, at2)po}- 

(41) 

This complete the derivation of Eq. (251. 

As a concrete example, let us work out the imple¬ 
mentation of Co explicitly. For this, we need to fix 
9 = 7r/4 in 6 = cos(0)ai -I- sin(0)a2 throughout the se¬ 
quence. It suffices to consider Tj = t for each Ram¬ 
sey measurement inducing the same required length for 
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the displacement amplitude. Therefore, in the first mea¬ 
surement ati = a(r)e“*T The desired amplitude is 
ati = V^a so that 'Db{at^) = Vai{a)'Da^[a.). The only 
thing we need to fix is |a| = 1.34. The second mea¬ 
surement occurs at time t 2 giving at^ = such 

that uj{t 2 — ti) = 27r/3. Therefore, we have 'Db{at 2 ) = 
/^). In the third measurement, 
at 3 = a(r)e“*3. Likewise, we chose w(t 3 — t 2 ) = 2tt/S 
so that Vbiat^) = 1'^'). With 

this, we could implement the required displacement op¬ 
erators appeared in C 3 using the sequence of Ramsey 
measurements. Generally, different three-point correla¬ 
tions of the associated modular variables composing ( 10 ) 
are obtained by suitable choice of ip \, (^2 and (^3 in each 
sequence of Ramsey measurements. Similar considera¬ 
tions are used for implementing other parts of the PM 
square. 


III. FORCED QUANTUM HARMONIC 
OSCILLATOR 


In this section we briefly summarize the mechanisms 
for creating large enough spin-dependent displacements 
of vibrational modes as it is needed for the contextual- 
ity tests proposed in this work. The basic interaction 


Hamiltonian as given in Eq. (22 1 can be realized in a 


wide variety of systems, ranging from trapped ions to 
nanomechanical resonators. For example, in the case of a 
trapped ion the hyperfine levels are coupled to the vibra¬ 
tional modes of the ion via state dependent optical dipole 
forces, which can be modulated at different frequencies to 
excite selective modes [SSUMl. Equivalent state depen¬ 
dent forces can be implemented by coupling the motion 
of a trapped ion or a nanomechanical resonator to a spin 
1/2 system via strong magnetic field gradients [T7U3^ . In 
this case a time-dependent coupling can be generated by 
simply flipping the spin state periodically. The common 
aspect is a modulated-in-time coupling strength which 
can resonantly drive the harmonic oscillator to an am¬ 
plitude which is much larger than what can be achieved 
with a static coupling. To illustrate this we consider the 
Hamiltonian 


= A(t)5l'-kA*(t)6, (42) 


ignoring the state-dependence of the force for now. The 
time evolution in the interaction picture generated by 
this Hamiltonian for duration t = t — Iq is obtained by 
solving 

dtU{t) = -i[A(t) 6 l'e*‘^‘ -p \*{t)be-^‘^^]U{t). (43) 

We characterize the solution by the Ansatz U{t) = 
with d(r) defined as 


d(r) 


-i [ A(<')e*‘^*'dt', 

•J tn 


(44) 




FIG. 4: Evolution of the displacement amplitudes a{r) (left) 
and q(' 7") (right) when excited with a near resonantly oscil¬ 
lating force. The upper panels correspond to the case where 
A(t) = The lower panels correspond to the case 

where A(t) = Ao cos[(a; -I- 5)t]. 


and 


(j){T) = C dt' r dt" sin(a;(t' - t")). (45) 

J tc\ J tn 


f to ^ to 

Therefore, the time evolution in Schrodinger picture is 


(46) 


Alternatively, we can rewrite (46) as 

[/(r) = 


(47) 


in which case a^r) = e“®‘^’’d(r). 

Let us give examples of two different types of cou¬ 
pling by which large enough displacement of a vibra¬ 
tional mode can be achieved. One example is X(t) = 
where oj is the natural frequency of an os¬ 
cillator and S = ujQ — uj characterizes the detuning of 
the driving (laser) field frequency from the vibrational 
frequency. Therefore, we have 

d(r) = -i f Aoe*'5*'dt' = ^(1 - e*^"), 

Jto ° (4g) 

<()(r) = sindt). 

Let us consider another type of interaction character¬ 
ized by a real coupling strength X{t) = Aocos[(a; -I- (5)t], 
where ui is the natural frequency and d = wq — w is the 






















detuning. We have 


t 

Ao cos[(a; + (5)t']e*“* dt'. (49) 

^0 

This finally (with defining a(T) = aa^r) + ia/(r)) yields 

= sis/t+u) 

and 

a/(^) = ^ 2a;) sin(r((5 + w)) - w sin(a;r)]. 

(51) 

The corresponding displacement amplitudes are shown in 

Fig-IH 


[1] E. P. Specker, Dialectica 14, 239 (1960). English trans¬ 
lation: arXiv:1103.4537. 

[2] S. Kochen, and E. P. Specker, J. Math. Mech. 17, 59 
(1967). 

[3] A. A. Klyachko, M. A. Can, S. Biniciogln, and A. S. 
Shumovsky, Phys. Rev. Lett. 101, 020403 (2008). 

[4] J. S. Bell, Physics (Long Island City, N.Y.) 1, 195 (1964). 

[5] A. Cabello, Phys. Rev. Lett. 101, 210401 (2008). 

[6] P. Badzi§,g, 1. Bengtsson, A. Cabello, and I. Pitowsky, 
Phys. Rev. Lett. 103, 050401 (2009). 

[7] S. Yu and C. H. Oh, Phys. Rev. Lett. 108, 030402 (2012). 

[8] M. Kleinmann, C. Budroni, J.-A. Larsson, O. Giihne, 
and A. Cabello, Phys. Rev. Lett. 109, 250402 (2012). 

[9] R. Ramanathan and P. Horodecki, Phys. Rev. Lett. 112, 
040404 (2014). 

[10] A. Cabello, M. Kleinmann, and C. Budroni, 
arXiv:1501.03432. 

[11] G. Kirchmair, F. Zahringer, R. Gerritsma, M. Klein¬ 
mann, O. Giihne, A. Cabello, R. Blatt, and C. F. Roos, 
Nature (London) 460, 494 (2009). 

[12] O. Moussa, C. A. Ryan, D. G. Cory, and R. Laflamme, 
Phys. Rev. Lett. 104, 160501 (2010). 

[13] E. Amselem, M. Radmark, M. Bourennane, and A. Ca¬ 
bello, Phys. Rev. Lett. 103, 160405 (2009). 

[14] J. Ahrens, E. Amselem, A. Cabello, and M. Bourennane, 
Sci. Rep. 3, 2170 (2013) 

[15] R. Lapkiewicz, P. Li, C. Schaeff, N. Langford, S. 
Ramelow, M. Wiesniak, and A. Zeilinger, Nature (Lon¬ 
don) 474, 490 (2011) 


[16] X. Zhang, M. Um, J. Zhang, S. An, Y. Wang, D. Deng, 
C. Shen, L.-M. Duan, and K. Kim, Phys. Rev. Lett. 110, 
070401 (2013). 

[17] M. Howard, J. J. Wallman, V. Veitch, and J. Emerson, 
Nature (London) 510, 351 (2014). 

[18] R. Raussendorf, Phys. Rev. A 88, 022322 (2013) 

[19] N. Delfosse, P. A. Guerin, J. Bian, and R. Raussendorf, 
arXiv: 1409.5170. 

[20] K. Horodecki, M. Horodecki, P. Horodecki, R. Horodecki, 
M. Pawlowski, and M. Bourennane, arXiv: 1006.0468. 

[21] O. Giihne, G. Budroni, A. Cabello, M. Kleinmann, and 
J.-A. Larsson, Phys. Rev. A 89, 062107 (2014). 

[22] G. Ganas, M. Arias, S. Etcheverry, E. Gomez, A. Cabello, 
G. B. Xavier, and G. Lima, Phys. Rev. Lett. 113, 090404 
(2014). 

]23] K. Banaszek and K. Wodkiewicz, Phys. Rev. A 58, 4345 
(1998) 

]24] A. R. Plastino and A. Cabello, Phys. Rev. A 82, 022114 

( 2010 ) 

]25] S. Massar and S. Pironio, Phys. Rev. A 64, 062108 

( 2001 ) 

[26] A. Asadian, C. Brukner, P. Rabl, Phys. Rev. Lett. 112, 
190402 (2014) 

[27] C. Gneiting and K. Hornberger, Phys. Rev. Lett. 106, 
210501 (2011). 

[28] G. Zhang, S. Yu, Q. Chen, C.H. Oh, Phys. Rev. Lett. 
Ill, 190501 (2013). 

]29] A. Peres, Phys. Lett. A 151, 107 (1990). 

[30] N.D. Mermin, Phys. Rev. Lett. 65, 3373 (1990) 

[31] O. Giihne, M. Kleinmann, A. Gabello, J.-A Larsson, G. 
Kirchmair, F. Zahringer, R. Gerritsma, and G. F. Roos, 
Phys. Rev. A 81, 022121 (2010) 

[32] J. Szangolies, M. Kleinmann, and O. Giihne, Phys. Rev. 
A 87, 050101(R) (2013). 

]33] A Vourdas, Rep. Prog. Phys. 67, 267 (2004) 

[34] S.-Z. Zhang, X.-T. Xie, and W.-X. Yang, Gommunica- 
tions in Theoretical Physics 46, 306 (2006). 

[35] G. Monroe, D. Meekhof, B. King, D. Wineland, Science 
272, 1131 (1996). 

]36] P. C. Haljan, K.-A. Brickman, L. Deslauriers, P. J. Lee, 
and G. Monroe, Phys. Rev. Lett. 94, 153602 (2005). 

[37] F. Mintert and C. Wunderlich, Phys. Rev. Lett. 87, 
257904 (2001). 

]38] K. Lake, S. Weidt, J. Randall, E. D. Standing, S. G. 
Webster, and W. K. Hensinger, Phys. Rev. A 91, 012319 
(2015). 

]39] P. Rabl, P. Cappellaro, M. V. Gurudev Dutt, L. Jiang, J. 
R. Maze, and M. D. Lukin, Phys. Rev. B 79, 041302(R) 
(2009). 








